Multichannel calculation of excited vector cj) resonances and the 0(2170) 
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A multichannel calculation of excited J PC = 1 states is carried out within a generalization 
of the Resonance-Spectrum Expansion, which may shed light on the classification of the 0(2170) 
resonance, discovered by BABAR and originally denoted X (2175). In this framework, a complete 
spectrum of bare ss states is coupled to those OZI-allowed decay channels that should be most 
relevant for the considered energy range. The included S- and P-wave two-meson channels comprise 
the lowest pseudoscalar, vector, scalar, and axial-vector mesons, while in the qq sector both the 3 Si 
and 3 Di states are coupled. The only two free parameters are tuned so as to reproduce mass and 
width of the 0(1020), but come out reasonably close to previously used values. Among the model's T- 
matrix poles, there are good candidates for observed resonances, as well other ones that should exist 
according to the quark model. Besides the expected resonances as unitarized confinement states, a 
dynamical resonance pole is found at (2186 — i246) MeV. The huge width makes its interpretation as 
the 0(2170) somewhat dubious, but further improvements of the model may change this conclusion. 

PACS numbers: 14.40.Cs, 11.80.Gw, 11.55.Ds, 13.75.Lb 



I. INTRODUCTION 

In 2006, the BABAR Collaboration announced [J 
the discovery of a new vector-meson resonance, called 
X(2175), in the initial-state-radiation process e + e~ — > 
K+K-jtttj, observed in the channel 0(1O2O)/ O (98O), 
with the meson decaying to K + K~ and the /o(980) to 
7r + 7r~ or 7T°7r . Two years later, the BES Collaboration 
confirmed [|J this resonance, then denoted y(2175), in 
the decay J/t/j -> ry[-^ 77 ] <f>[-+K+K-] / (980)[->7r+7r-]. 
At present, the new state is included in the PDG list- 
ings as the 0(2170) though not in the summary ta- 
ble, with average mass M = (2175 ± 15) MeV and width 
r = (61±18) MeV. However, these resonance parameters 
are being strongly challenged by the very recent Belle Q 
results on the F(2175), alias 0(2170), observed in the pro- 
cess e+e- 07T+7T-, yielding M = (2079 ± 13±^|) MeV 

and T = (192 ± 23^) MeV. 

The observation of this highly excited 0-type resonance 
with (probably) modest width, besides the peculiar, 
seemingly preferential, decay mode 0/ o (98O), triggered 
a variety of model explications, most of which proposing 
exotic solutions. Let us mention first a strangeonium- 
hybrid (ssg) assignment, in the flux-tube as well as the 
constituent-gluon model [|[, and a perturbative compar- 
ison of 0(2170) decays in these exotic ansatzes with a 
standard 2 3 Di ss description from both the flux-tube 
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and the 3 P model, by the same authors [||. Other ap- 
proaches in terms of exotics, with QCD sum rules, are 
an ssss tetraquark assignment 0, and an analysis Q 
exploring both ssss and ssss configurations. In an effec- 
tive description based on Resonance Chiral Perturbation 
Theory 0], the bulk of the experimental data is repro- 
duced except for the 0(2170) peak. This then led to a 
3-body Faddeev calculation [10| . with the pair interac- 
tions taken from the chiral unitary approach. Indeed, a 
resonance with parameters reasonably close to those of 
the 0(2170) is thus generated, though a little bit too nar- 
row. Finally, a review on several puzzling hadron states 
PH mentions the possibility that the 0(2170) arises from 
S'-wave threshold effects. 



In the present paper, we shall study the possibility that 
the 0(2170) is a normal excited meson, by coupling a 
complete confinement spectrum of ss states to a variety 
of S- and P-wave two-meson channels, composed of pairs 
of ground-state pseudoscalar (P), vector (V), scalar (S), 
and axial- vector (A) mesons. The employed formalism is 
a multichannel generalization of the Resonance-Spectrum 
Expansion (RSE) [HG3, 

which allows for an arbitrary 
number of confined and scattering channels [l4| . 



In Sec. 2 the RSE is very briefly reviewed and the ex- 
plicit T-matrix for the present coupled-channel system is 
given. Resonance poles and their trajectories are shown 
in Sec. 3, and model cross sections in Sec. 4. We draw our 
conclusions in Sec. 5 and discuss possible improvements. 
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II. THE RSE APPLIED TO RECURRENCES 

The RSE model has been developed for meson-meson 
scattering in non-exotic channels, whereby the interme- 
diate state is described via an infinite tower of s-channel 
qq states. For the spectrum of the latter, in princi- 
ple any confinement potential can be employed, but in 
practical applications, a harmonic oscillator (HO) with 
constant frequency has been used, with excellent re- 
sults. For more details and further references, see Refs. 

mm m Emm. 

In the present investigation of strangconium vector 
mesons, both the 3 Si and 3 Di ss confinement channels 
are included, to be compared with recent work re- 
stricted to the 3 S± component only. We could in principle 
also consider deviations from ideal mixing, by coupling 
the corresponding two (uu + dd)/V2 channels as well, 
but such fine corrections will be left for possible future 
studies. For the meson-meson channels, we consider the 
most relevant combinations of ground-state P, V, S, and 
A mesons that have nonvanishing coupling to either of 
the two confinement channels in accordance with the 3 P 
model and the OZI rule. The resulting 17 channels are 
listed, with all their relevant quantum numbers, in Ta- 
ble |TJ For the channels containing an rj or rf meson, we 
assume a pscudoscalar mixing angle of 37.3°, in the flavor 
basis, though our results are not very sensitive to the pre- 
cise value. Also note that channels with the same parti- 
cles but different relative orbital angular momentum L or 
total spin S are considered different. This is only strictly 
necessary for different L, because of the corresponding 
wave functions, but is also done when S is different, for 
the purpose of clarity. All relative couplings have been 
computed using the formalism of Ref. [20[ , based on over- 
laps of HO wave functions. They are given in Table Q] for 
the lowest recurrences (n = 0). As a matter of fact, we 
list their squares, which are rational numbers, but given 
as rounded floating-point numbers in the table, also for 
clarity's sake. For higher n values, the couplings fall off 
very rapidly. Their n dependence, for the various sets 
of decay channels, is presented in Table ITT1 The thresh- 
old values in Table U are obtained by takingthe meson 
masses given in the PDG tables or listings with the 
exception of the K q (800) (alias k), for which we choose 
the real part of the pole position from Ref. [16| . as it 
lies closer to the world average of k masses. Note that 
we take sharp thresholds, even when (broad) resonances 
arc involved. We shall come back to this point in the 
conclusions. Finally, we should notice that a number of 
channels that also couple to ss vector states according to 
the scheme of Ref. [20J, viz. P-wave channels involving 
axial- vector mesons as well as some channels with tensor 
mesons, have not been included in the final calculations 
presented here. However, their influence has been tested 
and turned out to be very modest, due to the correspond- 
ing small couplings. 

Coming now to the explicit expressions for our model, 
let us first write down the effective meson-meson interac- 



tion, which consists of an intermediate-state s-channel qq 
propagator between two quark-antiquark-meson-meson 
vertex functions for the initial and final state, reading 

HQ 

V^ipuP^E) = A 2 /; ■ l ,,u;R. ll U-Aji :,/„; , (1) 
with 

^ae)= E E ) . ( 2 ) 

Z c =0,2n=0 & — £j n 

where the RSE propagator contains an infinite tower of 
s-channel bare qq states, corresponding to the spectrum 
of an, in principle, arbitrary confining potential. Here, 
En ^ is the discrete energy of the n-th recurrence in the 
ss channel with angular momentum l c , and n ^ is the 
corresponding coupling to the i-th meson-meson chan- 
nel. Furthermore, in Eq. ([T]), A is an overall coupling, 
and ji,.(pi) and pi are the L^-th order spherical Bessel 
function and the (relativistically defined) off-shell rela- 
tive momentum in meson- meson channel i, respectively. 
The spherical Bessel function originates in our string- 
breaking picture of OZI-allowcd decay, being just the 
Fourier transform of a spherical delta function of radius 
a. Together with the overall coupling constant A, the ra- 
dius a is a freely adjustable parameter here, though its 
range of allowed values turns out to be quite limited in 
practice. The couplings g^ n -. in Eq. @ arc obtained by 
multiplying the ones in Table U by those in Table UU for 
the corresponding channels. Because of the fast decrease 
of the latter for increasing n, practical convergence of the 
infinite sum in Eq. @ is achieved by truncating it after 
20 terms. 

Because of the separable form of the effective meson- 
meson interaction in Eq. (TT]), the fully off-shell T-matrix 
can be solved in closed form with straightforward algebra, 
resulting in the expression 

T^'^ipi^E) = -2aA 2 V^jl i (p 4 a)x 

N 

£ { t 1 - n n ^ } mj 4 (pJ°) > ( 3 ) 

m— 1 

with 

Qij(kj) = —2ia\ 2 fijkjj 3 Lr (kja)h^ : '(kja)Sij , (4) 

where (kja) is the spherical Hankel function of the 
first kind, kj and fj,j are the on-shell relative momen- 
tum and reduced mass in meson-meson channel j, re- 
spectively, and the matrix H(E) is given by Eq. ([2"|). 

As mentioned above, we assume an HO confinement 
spectrum with constant frequency, given by 

E n = 2m s + u)(2n + l c + 1.5) , (5) 
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TABLE I: Included two-meson channels, their internal and 
relative angular momenta and spins, couplings squared for 
n = 0, and thresholds. See Ref. 0) for properties of listed 
mesons, except for the Kg (800), discussed in the text. 



Channel 


9fl a =0,n) x 4" 


5a„=2,«) x 4" 


PP 

PV 
VV 

sv 

PA 
VA 


(2n + 3)/3 
(2n + 3)/3 
(2n + 3)/3 
(2n-3) 2 /9 
(2n - 3) 2 /9 
(2n-3) 2 /9 


71+1 

n + 1 

n + 1 
(n+l)(2n + 5)/5 
(n + l)(2n + 5)/5 
(n + 1) 2 



TABLE II: Dependence of couplings squared on recurrence n. 



with w = 190 MeV and m s = 508 MeV fixed at values 
determined long ago [2l[ . Some of the resulting bare HO 
ss states are given in Table Mil 

Now that the model has been fully defined, we are in 
a position to evaluate the on-shell components of the T- 
matrix defined in Eqs. (|3l2p . for the channels given in 
Tables HHnU 



III. EXPERIMENTAL STATUS OF STATES 

Before adjusting our two free parameters A and a from 
Eq. ([3]), let us first have a look at the experimental status 
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1301 
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2441 


2821 



TABLE III: Masses of bare ss states in MeV, for HO potential 
with w = 190 MeV and m s = 508 MeV (see Eq. (f5j and 
Ref. HI). 





M (MeV) 


T (MeV) 


0(1020) 
0(1680) 
0(2170) 


1019.455±0.020 
1680±20 
2175 ± 15 


4.26±0.04 
150±50 
61 ± 18 



TABLE IV: Listed J PC = 1~ resonances, with masses and 
widths [| (values for 0(1680) are estimates H). 



of vector resonances. According to the 2008 PDG list- 
ings 3, there are only 3 observed states, viz. the </)(1020), 
0(1680), and 0(2170), with the latter resonance omitted 
from the summary table. Their PDG masses and widths 
are given in Table IIVI Clearly, this is a very poor sta- 
tus, as several additional states must exist in the energy 
range 1-2 GeV according to the quark model, and also 
if we compare with e.g. observed p resonances Q in the 
same energy interval. Moreover, the 0(1680) can hardly 
be the first radial excitation of the 0(1020), in view of 
the well established A"*(1410), which is almost 300 MeV 
lighter, and a typical mass difference of 100-150 MeV 
between the strange and nonstrange (u, d) constituent 
quarks [U [22j ■ This conclusion is further supported if 
indeed the ,0(1250) is confirmed as the first radial recur- 
rence of the p(770) HHH. So the 0(1680) is more likely 
to be the 1 3 Di state, with a hitherto undetected 2 3 5i 
state somehwere in the mass range 1.5-1.6 GeV. As a 
matter of fact, in Ref. [24| a vector resonance was re- 
ported at roughly 1.5 GeV, though this observation is, 
surprisingly, included under the 0(1680) entry Even 
more oddly, another 0-like state, at ~ 1.9 GeV and re- 
ported in the same paper [24|, is also included under 
the 0(1680) @. However, a resonance at about 1.9 GeV 
should be a good candidate for the next radial ss recur- 
rence, if we take the observed p resonances in Ref. [23| 
for granted. 



IV. HUNTING AFTER POLES 

In view of the poor status of excited states, let us 
adjust our parameters A and a to the mass and width 
of the 0(1020). Here, we should mention that an ad- 
ditional phenomenological ingredient of our model is an 
extra suppression of subthreshold contributions, using a 
form factor, on top of the natural damping due to the 
spherical Bcssel and Hankcl functions in Eq. Such 
a procedure is common practice in multichannel phase- 
shift analyses. Thus, for closed meson- meson channels 
we make the substitution 



(ff(/c,n)) (S(j ,n)) 



for 3?efc?<0 



(6) 



The parameter a is chosen at exactly the same value as 
in previous work [l6|, [l8|, viz. a — 4 GeV -2 . 

Choosing now A = 4 GcV~ 3/2 and a = 4 GcV~\ 
we manage to reproduce mass and width of the 0(1020) 
with remarkable accuracy, namely = 1019.5 MeV 
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6 


2397 


-214 
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7 


2415 


-6 


2371 


-29 


conf., n = 3, 4 3 5"i 


8 






2415 


-8 


conf., n = 2, 3 3 Z>i 


9 


2501 


-236 


2551 


-193 


continuum 



TABLE V: Complex-energy poles in MeV, for 3 iSi ss channel 
only, and for both 3 5i and 3 Di. See text for further details. 



and = 4.4 MeV. Note that these values of A and a are 
of the same order of magnitude as in the work mentioned 
before [HI, HH, which dealt with scalar mesons. 

In Table IVl we collect all resonance poles encountered 
on the respective physical Ricmann sheets, which corre- 
spond to 3m ki > for closed channels and 3m ki < 
for open ones. When the latter conditions are not ful- 
filled, we call the corresponding Ricmann sheets unphys- 
ical. Moreover, we also show here the pole positions ob- 
tained by taking only the 3 S± ss channel and switching 
off the 3 D%, for fixed A and a. Focusing for the moment 
on those poles that originate in the states of the confine- 
ment spectrum (indicated by "conf." in the table), we 
see good candidates for the resonances at ~ 1.5 GeV and 
~ 1.9 GeV reported in Ref. (24|, and possibly also for 
the 0(1680), though our 1 3 D\ state seems somewhat too 
light. Note, however, that under the 0(1680) entry Q in 
the PDG listings there is a relatively recent observation 
[25l ] with a mass of (1623 ±20) MeV, which is compatible 
with our pole at 1602 MeV. Furthermore, the imaginary 
parts of the confinement poles are generally too small, 
except for the 0(1020). We shall come back to this point 
in the conclusions below. Besides the latter poles, also 
two so-called continuum poles are found, often designated 
as dynamical poles, the most conspicuous of which is the 
one at (2186 — i246) MeV, as the real part is very close to 
the mass of the 0(2170) as measured by BABAR [l| and 
BES @- However, in view of the much too large width, 
even as compared to the Belle Q value, considerable cau- 
tion is urged. Also this point will be further discussed in 
the conclusions. 

Some words are in place here about our identification 
of the 3 Si and 3 Di confinement poles in Table [V] The 
point is that, rigorously speaking, these designations only 
make sense for pure confinement states and, moreover, 
without any 3 Si / 3 D\ mixing. Now, in our approach, the 
very mixing is provided by the coupling to common de- 
cay channels. So for any nonvanishing value of the overall 
coupling A there are no longer pure 3 S\ and 3 Di states, 
while for the physical value of A the mixing is proba- 
bly considerable. Moreover, there is no obvious way to 
tell which pole of a pair originating in a degenerate con- 




Re E (GeV) 



FIG. 1: Trajectory of first continuum pole, for 2.26 < A < 
5.99 (GeV" 3/2 ), from left to right. Bullet represents A = 
4 GeV~ 3//2 , while dashed line indicates unphysical Riemann 
sheet. 

finement state stems from either 3 S\ or 3 D\. Therefore, 
our identification is partly based on the couplings in Ta- 
ble HI which on the whole suggest larger shifts for 3 S\ 
than for 3 D\ , partly on a comparison with a perturbative 
approach employed in Ref. [26| to find poles for small A. 

The designation continuum pole becomes clear when 
plotting a corresponding trajectory as a function of the 
overall coupling A. In Fig. [TJ the first such pole is 
shown to have an increasingly large imaginary part for 
decreasing A, eventually disappearing in the continuum 
for A — ► 0. Note that the small jump at the important 
S'-wave K*Ki(\270) threshold is due to a minor thresh- 
old discontinuity of the damping function in Eq. ([6|) for 
complex momenta. 

Figure [5] shows a similar trajectory, but now for the 
lowest confinement pole, which ends up as the 0(1020) 
resonance. Notice the large negative mass shift (« 
280 MeV), as well as the way the pole approaches the 
KK threshold, which is typical for P-wave decay chan- 
nels. Also note that the tiny jump in the trajectory is 
due to the way relativistic reduced mass is defined be- 
low threshold, which in the case of closed channels with 
highly unequal masses (KKi(Y27§) here) requires an in- 
tervention to prevent the reduced mass from becoming 
negative. 

In Fig. [31 we depict the trajectories of the 2 3 S\ and 
1 3 Di confinement poles. Note that the coupling to decay 
channels lifts the original degeneracy of the 2 3 S± and 
1 3 Di HO states. 

The trajectories of the next pair of confinement poles, 
i.e., 3 3 5i and 2 3 Di, are drawn in Fig.g] Note the highly 
nonlinear behavior of the poles, showing the unreliabil- 
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FIG. 2: l 3 ^! confinement pole for 4.31 > A > (GeV" 3/2 ) 
Bullet represents A = 4 GeV~ 3/2 . 




1.95 



2.00 

RcP (GeV) 



2.05 



FIG. 4: 3 3 Si (lower) and 2 3 Di (upper) confinement poles 
for 4.2 > A > (GeV _3/2 ) and 5.99 > A > (GeV _3/2 ), 
respectively. Bullets represent A = 4 GeV _3//2 , while dotted 
and dashed lines indicate unphysical Riemann sheets. 



ity of perturbative methods to estimate coupled-channel 
effects. 
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FIG. 3: 2 3 Si (lower) and 1 A D\ (upper) confinement poles 
for 5.0 > A > (GeV~ 3/2 ) and 4.76 > A > (GeV~ 3/2 ), 
respectively. Bullets represent A = 4 GeV~ 3,/2 , while dotted 
and dashed lines indicate unphysical Riemann sheets. 



1.2 1.4 

E (GeV) 

FIG. 5: Elastic P-wave KK cross section. Full line: both 3 5i 
and 3 Di ss channels included; dashed line: only 3 Si. 



V. CROSS SECTIONS 

Now we shall show, as mere illustrations, some of the 
cross sections related to the resonance poles found in the 
preceding section. In Fig. the clastic P-wave KK 
cross section is depicted in the energy region covering 
the ^(1020) as well as the 2 3 Si and 1 3 Di resonances. 
We see that including the 3 Di ss channel has the effect of 



6 




— i — i — i — i — i—i — i — i — i — | — i — i — i — i — | — i — i — i — i — i 
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FIG. 6: Natural logarithm of the ratio of the elastic KK* and 
KK cross sections. 



lowering the 2 3 Si state, besides the generation of an ad- 
ditional resonance, of course. This "repulsion" between 
the 3 5i and 3 D\ poles is also noticed for the 3 3 Si and 
2 3 Di states. 

Figure [5] shows the relative importance of the KK and 
KK* channels in the energy interval 1.5-1.7 GeV, which 
should be relevant for the 0(1680). The plotted quan- 
tity is the logarithm of the ratio of the elastic KK* and 
KK cross sections, which shows that the KK* channel 
is strongly dominant, except at low energies, because of 
phase space, and close to the pole at ~ 1.6 GeV, where 
the two cross sections are comparable. Dominance of the 
KK* decay mode is reported under the 0(1680) PDG 
entry 

Turning now to the 0(2170) energy region, we show 
in Fig. [3 the elastic S- and D-wave 0(1O2O)/ O (98O) 
cross sections. The effect of the continuum pole at 
(2186 — z246) MeV is noticeable as a small and very broad 
enhancement in the D-wave cross section. In the S-wave 
case, its effect is completely overwhelmed by the huge 
cross section at threshold, partly due to the 3 3 Si pole 
not far below. Also quite conspicuous are the here pre- 
dicted 4 3 iSi and 3 3 D\ resonances (see Table ED for the 
respective pole positions). Of course, all these model 
elastic cross sections have little direct bearing upon the 
experimentally observed production cross sections. The 
production process of the 0(2170) may be studied with 
the RSE production formalism [23], but that lies out- 
side the scope of the present investigation, which focused 
on the possibility of generating a 0(2170) resonance pole 
through coupled channels. 

Finally, in Fig. [8] we plot the logarithm of the ratios 
of the elastic S-wave 0(lO2O)/o(98O) cross section and 
the elastic K*K*, 0(1020)?/, and K*K X {1270) cross sec- 




2.1 2.2 2.3 2.4 

E (GeV) 



FIG. 7: Elastic D-wave (solid line) and S'-wave (dashed line) 
0(1020) fo (980) cross section. Dotted line: S'-wave cross sec- 
tion for 3 Si channel only. 



tions, in the energy interval 2.0-2.3 GeV. We see that 
the S'-wave 0(lO2O)/o(98O) cross section dominates up 
to about 2.08 GeV, but getting overwhelmed first by the 
(P-wave) K*K* channel, and then even more so by the 
S'-wave K* Ki(1270) channel, right from its threshold at 
«2.16 GeV upwards. Also the 0(lO2O)7y' channel is be- 
coming more important here. As for the K*K* channel, 
it gives rise to a final state with two kaons and two pions, 
i.e., the same as that for which the 0(2170) was observed. 
So the experimental status of the 0(2170) might be im- 
proved if one succeeded in identifying and isolating the 
K*[— ► Ktt]K*[-^ Ktt] decay mode, which should be quite 
important. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we have applied the RSE formalism for 
non-exotic multichannel meson-meson scattering to cal- 
culate the resonance spectrum of excited vector mesons, 
and to find out whether this way the 0(2170) can be gen- 
erated. The inclusion of all relevant two-meson channels 
that couple to the bare 3 Si and 3 Di ss states should 
guarantee a reasonable description. Thus, several vector 
resonances are predicted, some of which are good can- 
didates for observed states, while others may correspond 
to others, undetected so far, but quite plausible in view 
of observed partner states in the excited p spectrum. Fi- 
nally, a very broad 0-like resonance pole of a dynamical 
origin is found, with real part very close to that of the 
0(2170), but a much too large imaginary part, so that 
its interpretation remains uncertain. On the other hand, 
the calculated resonances originating in the confinement 
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2.10 2.20 
E (GeV) 

FIG. 8: Natural logarithm of the ratios of the elastic S- 
wave <f>{ 1020) fo (980) cross section and the elastic K*K* (solid 
line), 0(lO2O)r/ (dotted line), and (1270) (dashed line) 

cross sections. 

spectrum are generally too narrow. 

These considerations bring us to the main problem of 
our description, namely the inclusion of sharp thresh- 
olds only. The point is that many of the channels in Ta- 
ble U involve highly unstable particles, several of which 
are broad to very broad resonances themselves. Treating 
the corresponding thresholds as sharp is clearly an ap- 
proximation. In particular, the /o(980) meson included 
in the 0(1020) fo (980) channels is a very pronounced res- 
onance in the coupled tttt-KK system. This feature is 
crucial in the three-body calculation of the 0(2170) in 
Ref . [l(| , which indeed produces a clear resonance signal 
at almost the right energy, and even with a somewhat 



too small width. We believe that in our approach, too, a 
narrower 0(2170) might be generated, if we could account 
for the physical width of the /o(980) meson, and also for 
the widths of the K* and Xi(1270) resonances in the 
here included if*ifi(1270) channel. The reason is that 
the widths effectively cause these channels to act already 
below their central thresholds, which will strongly influ- 
ence poles just underneath. Especially the width of very 
strongly coupling K* Ki(l270) channel, whose threshold 
lies only some 25 MeV below the real part of the con- 
tinuum pole at (2186 — i x 246) MeV, will surely have a 
very significant effect on this pole's trajectory. Because of 
the typical behavior of continuum poles, with decreasing 
width for increasing coupling, we expect that the width of 
our 0(2170) candidate may thus be reduced. Conversely, 
including the widths of final-state resonances will proba- 
bly increase the widths of the now too narrow excited 
resonances stemming from the confinement spectrum. 

Of course, to account for the nonvanishing widths 
of mesons in the coupled channels is a very difficult 
problem, since the simple substitution of the here used 
real masses by the true complex masses will destroy the 
manifest unitarity of the S'-matrix. Work is in progress 
to redefine the S'-matrix for such cases, so as to enforce 
its unitarity by construction. 
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